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Summary 

Nonlinear finite-element computer programs are too 
costly to use in the early design stages for hot-section 
components of aircraft gas-turbine engines. To improve 
the durability of these components, it is necessary to 
develop simpler and more economical methods for 
representing the structural response of materials under 
cyclic loading. Programs have been under way at the 
NASA Lewis Research Center to develop simplified 
procedures for performing nonlinear structural analysis 
using only an elastic finite-element solution or strain- 
gauge data as input. 

Development of the simplified method was based on 
the assumption that the inelastic regions in the structure 
are constrained by the surrounding elastic material. This 
implies that the total strain history can be defined by an 
elastic analysis. Initial development of the method did 
not account for any strain redistribution. A computer 
program (ANSYMP) was created to predict the stress- 
strain history at the critical fatigue location of a 
thermomechanically cycled structure from elastic input 
data. Appropriate material, stress-strain and creep 
properties, and plasticity hardening models were 
incorporated into the program. Effective stresses and 
equivalent plastic strains are approximated by an iterative 
and incremental solution procedure. Creep effects can be 
obtained on the basis of stress relaxation at constant 
strain, cumulative creep at constant stress, or a 
combination of stress relaxation and creep accumulation. 

The simplified procedure predicts the critical location 
of stress-strain response with reasonable accuracy relative 
to nonlinear finite-element analyses for thermally cycled 
problems. However, the limitation of no strain 
redistribution is most likely to be violated in the case of 
mechanically loaded structures, especially in the vicinity 
of stress concentrations. Nonlinear finite-element 
analyses of notched-plate specimens subjected to cyclic 
mechanical loading have shown that the total strain range 
at the critical location can be significantly larger than 
would be predicted from elastic solutions. 

This study derived and incorporated corrections in the 
simplified procedure to account for local total strain 
redistribution and residual stresses due to mechanical 
load cycling. These corrections would remove some of 
the limitations and extend the applicability of the 


procedure. The corrections were based on the Neuber 
rule relating the theoretical stress concentration factor to 
the actual stress and strain concentration factors. 

Two variations of a benchmark notched-plate problem 
were analytically examined to verify the accuracy of the 
improved simplified method. Verification was made 
through comparison with three-dimensional nonlinear, 
finite-element analyses. Cyclic stress-strain and creep 
properties for Inconel 718 alloy, a kinematic hardening 
model, and the von Mises yield criterion were used for the 
benchmark notched-specimen problems. Elastic and 
elastic-plastic finite-element analyses were performed by 
using the MARC nonlinear finite-element computer 
code. The elastic solutions for the critical locations were 
used as input data for the simplified analysis computer 
code. The stress-strain histories at the critical locations 
from the simplified and nonlinear finite-element analyses 
were compared. 

The comparisons demonstrated that the improved 
simplified method can duplicate the cyclic stress-strain 
hysteresis loops from the MARC elastic-plastic analyses 
to a high degree of accuracy. For the benchmark 
problem, ANSYMP used 0.3 percent of the central 
processor unit time required by MARC to compute the 
inelastic solution. 

Introduction 

The severe operating conditions in advanced aircraft 
gas-turbine engines have subjected hot-section 
components such as turbine blades and combustor liners 
to thermomechanical load cycles that induce repeated 
inelastic strains and eventual fatigue cracking. Life 
prediction and improvements in the durability of these 
components require accurate knowledge of the 
temperature stress-strain history at the critical crack 
initiation location of the structure. 

Nonlinear finite-element computer codes are 
increasingly being used for calculating inelastic structural 
response. However, nonlinear finite-element analysis is 
not generally feasible for use as a component design tool 
because of the high computing costs associated with the 
iterative and incremental nature of the inelastic solutions. 
Computing costs are further increased by the presence of 
high thermal gradients and geometric irregularities, such 
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as cooling holes. Most hot-section components are of 
such geometric complexity as to necessitate three- 
dimensional analyses, frequently with substructuring. 
Three-dimensional, nonlinear finite-element analyses are 
prohibitively time consuming and expensive to conduct in 
the early design stages for combustor and turbine 
structures. To improve the design of hot-section 
components, it is necessary to develop simpler and more 
economical methods for representing structural behavior 
under cyclic loading. 

The initial development of a simplified analytical 
procedure for estimating the stress-strain history of a 
thermomechanieally loaded structure subject to cyclic 
inelasticity was implemented in a computer program 
(ANSYMP) (refs. 1 to 3). This procedure predicts the 
stress-strain history at the critical location for crack 
initiation by using as input the total strain history 
calculated from elastic finite-element analyses. An 
incremental and iterative procedure estimates the plastic 
strains from the material stress-strain properties and a 
plasticity hardening model. Creep is incorporated in the 
program through options which permit the calculation of 
stress relaxation at constant strain, creep at constant 
stress, or a combination of stress relaxation and creep 
accumulation. 

Analytical predictions from the simplified method for 
a number of problems were compared in references 1 to 3 
with nonlinear finite-element solutions from the MARC 
computer program (ref. 4). These problems involved 
uniaxial- and multiaxial-stress states, isothermal and 
nonisothermal conditions, and various materials and 
plasticity hardening models. The problems included an 
Inconel 718 benchmark notched specimen that was 
mechanically load cycled in an experiment to verify 
structural analysis methods (ref. 5). In these analyses the 
predicted stress-strain response from the simplified 
procedure was in good agreement with the MARC 
solution for the thermally cycled problems but was 
truncated for the benchmark notch problem. This 
disparity in the latter case was due to strain redistribution 
under mechanical load cycling which resulted in a much 
larger total strain range than was predicted from the 
elastic finite-element analysis. 

In the present study, the simplified method was further 
developed to consider strain redistribution effects from 
plastic strain reversal under mechanical load cycling. 
Corrections were derived to account for the 
redistribution of local total strains due to plastic yielding 
and residual strains during load reversal. These 
corrections were incorporated into a version of the 
ANSYMP elastic-plastic-creep program. The improved 
procedure was applied to the benchmark notch problem. 
Cyclic stress-strain properties for Inconel 718 alloy, a 


kinematic hardening model, and the von Mises yield 
criterion were used for these analyses. Elastic and elastic- 
plastic finite-element analyses were performed by using 
the MARC nonlinear finite-element computer code. The 
elastic solution at the critical location was used as the 
total strain history input for the simplified analysis 
computer code. The improved simplified procedure was 
verified on the basis of how well it was able to duplicate 
the stress-strain hysteresis loops from the MARC elastic- 
plastic analyses of the benchmark notch problem. 

Symbols 

E modulus of elasticity 

E e modified modulus of elasticity 

E p work hardening slope (fig. 1) 

e nominal total strain at net section 
K,n temperature-dependent constants in cyclic 
stress-strain equation 

K ( theoretical elastic stress or strain concentration 

factor 

K e actual strain concentration factor 

K a actual stress concentration factor 

S effective nominal elastic stress at nth increment 
S* m maximum effective local elastic stress 

S* effective local elastic stress (real or imaginary) 
e equivalent total strain, unmodified 

e e uncorrected equivalent elastic strain 

€p equivalent plastic strain 

maximum equivalent plastic strain in cycle 
(fig- 1) 

e/? equivalent local residual strain 

e * equivalent total strain, modified 

fi Poisson’s ratio 

o effective stress (loading), unmodified 

o' maximum effective stress in cycle (fig. 1) 
o m maximum effective stress, unmodified 

o* m maximum effective stress, modified 

o R local residual stress 

o vi initial yield stress in loading part of cycle 

o'yj initial yield stress in unloading part of cycle 

02 effective stress (unloading), unmodified 

<72 effective stress (unloading), modilied 

o* effective stress (loading), modified 


2 


Analytical Procedure 


^Stress, o 


Improvements were made to a simplified inelastic 
procedure for calculating the stress-strain history at the 
critical fatigue location of a structure subjected to cyclic 
thermomechanical loading. Details of the simplified 
analysis procedure, including a computer program flow 
chart, are presented in reference 2. The basic assumption 
in the initial development of the procedure was that the 
inelastic region is local and constrained from 
redistribution by the surrounding elastic material. It 
follows from this assumption that the total strain history 
at the critical location can be defined by an elastic 
solution. Justification for the assumption of elastic 
constraint of local inelasticity can be found in references 
6 to 12, which present finite-element analyses of 
thermally cycled structures, such as combustor liners, 
aircooled turbine blades, and wedge-fatigue specimens, 
where the total strain ranges from elastic and nonlinear 
solutions are in close agreement. However, this 
assumption is invalid for structures subjected mainly to 
mechanical load cycling such as the benchmark notch 
problem (ref. 1). 

In this study, a strain redistribution factor (SRF) was 
derived to account for local total strain redistribution 
under applied cyclic loads. This correction also 
accounted for residual strains induced by plastic yielding. 
The strain redistribution factor is applied to the ideal 
local total strain obtained from the elastic solution. These 
corrections were incorporated in a version of the 
ANSYMP program by using a kinematic hardening 
model to characterize the yield surface under cycling. A 
representation of a cyclic stress-strain curve by a bilinear 
kinematic hardening model is illustrated in figure 1. The 
loci of the tips of the cyclic curves are described by the 
relation 
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The work hardening slope E p for the kinematic hardening 
model was determined from energy considerations to give 
the same strain energy, as indicated by the enclosed area 
in figure 1, as the actual stress-strain curve. This work 
hardening slope is defined by 



o' 2 n 
tp \ + n 


( 2 ) 


The von Mises criterion was used to convert the total 
strain from a uniaxial stress-strain curve into a modified 



equivalent total strain curve as discussed in reference 13. 
This modified elastic part of the equivalent total strain 
corresponds to the measured elastic strain multiplied by 
2(1 + ^)/3. This relationship must be taken into account 
in a multiaxial stress-strain field when applying results 
from elastic finite-element analyses. As a result, all 
stresses and strains from the elastic solution are expressed 
in terms of von Mises effective stress and modified 
equivalent total strain with signs assigned based on the 
dominant principal stress or strain direction. For 
computational convenience, the factor 2(l+/x)/3 is 
included in the elastic modulus, and the modified elastic 
modulus is designated by E e , where 


Ee ~2(\ +l x) E 


The elastic input data are subdivided into a sufficient 
number of increments to define the stress-strain cycle. 
Dwell times are specified for increments which require 
creep analysis. The increments are analyzed sequentially 
to obtain the cumulative plastic and creep strains and to 
track the yield surface. An iterative process is used to 
calculate the yield stresses for increments undergoing 
plastic straining. First, an estimated plastic strain is 
assumed for calculating an initial yield stress from the 
stress-strain properties and the simulated hardening 
model. Second, a new plastic strain is calculated as the 
difference between the total strain and the elastic- and 
creep-strain components. The yield stress is then 
recalculated by using the new plastic strain. This iterative 
process is repeated until the new and previous plastic 
strains agree within a tolerance of 1 percent. 

Neuber’s rule (ref. 14), which has been widely used to 
approximate local stresses and strains in the plastic 
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region, was selected as the basis for the development of 
the strain redistribution corrections. With minor 
modifications, this rule provides acceptable predictions 
of local stresses and strains in the loading part of a stable 
cycle and less accurate predictions in the unloading part 
of the cycle. The relative failure of Neuber’s rule for the 
unloading part of the cycle has been mainly due to neglect 
of the local residual stresses. Regardless of the part of the 
cycle being studied, an adjustment must be made to the 
stresses to account for the work hardening slope of the 
cyclic stress-plastic-strain curve. 

The strain redistribution analysis for the stable cycle is 
divided into the following three stages: elastic-plastic 
loading, elastic unloading, and inelastic unloading. 
Figure 2 is presented to further define the symbols used in 
the governing equations for each of these stages. 

Elastic-Plastic Loading 

Neuber’s rule specifies a relationship among the 
theoretical stress concentration factor K h the actual stress 
concentration factor K a , and the actual strain 
concentration factor K ( of the form 


A = 2 E p 

C = 4^5*2 

E e 

The local stresses and strains calculated by using the 
work hardening slope (eqs. (6) and (7)) have to be 
modified to account for the elastic-strain change between 
o and o V j. The modified local stress a* corresponding to 
the local strain c* can be expressed as 



and the modified local strain e* as 



K}=K a K ( 


The plastic strain e p at the stress level <x* is 
(3) 


Substituting local and net section stresses and strains into 
equation (3) gives 



(4) 


E e 


( 10 ) 


To relate the ideal local inelastic strain to the local total 
strain, a strain redistribution factor (SRF) is introduced. 
The SRF at the nth increment on the stress-strain curve 
can be written as 


where K t also equals S*/S or a = S* 2 /eE e . 

The cyclic stress-strain relation is assumed to take the 
form 


(SRF) n = e* 


Sn 

E e 


(ID 


a = E e t when (5) 

t e 

or 

° = ^ 1 - ^ ^ Oyi + Ep< when e > jjr (6) 

Combining equations (4) to (6) gives the local strain as 


To improve the accuracy of the predicted local stress- 
strain history, an SRF was generated for each increment 
during the loading part of the cycle. This SRF is the 
cumulative total strain correction up to the nth point on 
the stress-strain curve. Therefore, the total strain at that 
increment is the uncorrected elastic strain e e for that 
increment plus the sum of all the incremental SRF’s up to 
that increment: 


- B 4 - VB- + C 


(7) 


6m — 6mm 


£ [<SRF)j-(SRF)j_,l 

/=! L J 


( 12 ) 


Equation (12) is another way of expressing t * given in 
equation (9). 
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Elastic Unloading 

The local stress at the nth increment in the elastic 
unloading part of the cycle (fig. 2) is simply 


G n ~ G m ^n) 

and the corresponding local total strain is 


(13) 


£ n £ m 


S' m -S‘ n 


(14) 


or expressed in a form similar to equation (12) 


t j^(SRF); - (SRF)j _ i J 


(15) 


For an isothermal problem where an incremental 
approach is not required, equation (15) can be simplified 
to € n = 4- (SRF) m . 

Inelastic Unloading 

Upon complete removal of the applied load, residual 
stresses remain across the net section. Denoting the local 
residual stress as o R gives 


° R ~ °tn Sm ( 1 

The complimentary local residual strain, e R , is then 

( 17 ) 


°R 


e R = 


Again the cyclic stress-strain relation in the plastic 
region is assumed to take a linear form 




08 ) 


As a kinematic hardening model is used in the analysis, 
the new initial yield stress is given by 


°yi ~ 2 a v i 


(19) 


Considering the zero stress and plastic strain at the new 
origin o* m (fig. 2), the local stress and total strain can be 
expressed as 


<j-K 0 S+or 

and 

e=K ( e + t R 


( 20 ) 


( 21 ) 


Combining equations (3), (18), (20), and (21) leads to 
the solution for 


-A - V/42 + fl 

2E n 


(22) 


where 


' 4 =^ 1 - ^r^o'yi- EptR- °R 

[("ft)* 




5*2 

e 




} 


The actual local stresses and strains are obtained by a 
series of iterations between the modified and unmodified 
values. As shown here, two iterations are sufficient to 
give an acceptably convergent solution. Iterations are not 
necessary for the loading part of the cycle. The modified 
local stress corresponding to e* is expressed as 


a 5 = i> + ( 1 _ t ) <72 

where o 2 is obtained from the following steps: 


(23) 


It 

1 

°yi F Ep( 

(24) 

■*=&"W+| 

b e 

('-»• 

(25) 

a — a* 


(26) 

\~e — 



^ O yj + Epi 1 

(27) 
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and the modified local total strain e* is 


(28) 

(29) 

(30) 


* (7 — O'? O') — (7 5 

e * = e +e n 


(31) 


The plastic strain e p at the stress level is 


e = e* — — 
£p 6 E„ 


(32) 


The strain redistribution factor at the nth increment on 
the stress-strain curve is 


(SRF) n = e„ - 


°n ~ a n2 °n2 ~ a n2 $h + a R 


E e E e 


(33) 


or, in terms of the uncorrected elastic strain, the total 
strain at the nth increment can be written as 


e n - ten + (SRF) n 


n r 


(SRF)j-(SRF)j_ \ 


(34) 


where 


(SRF) m = e m — y? - - 


Equation (34) is simply another form of expressing the 
e* previously given in equation (31). 

Another version of the ANSYMP computer program 
was created to implement the improved simplified 
analytical procedure. The program consists of the main 
executive routine, ANSYMP, and four subroutines, 
ELAS, YIELD, CREEP, and SHIFT. The incremental 
elastic data and temperatures are read into subroutine 
ELAS. Material stress-strain properties as a function of 
temperature and a simulated hardening model are 
incorporated in subroutine YIELD, and the creep 


characteristics are incorporated in subroutine CREEP. 
Subroutine SHIFT is required to update the temperature 
effects on the yield stress shift. SHIFT also serves the 
function of deciding the future direction of the yield 
surface under nonisothermal conditions by determining 
the relation of future to past thermal loading. 

The calculation scheme initially follows the effective 
stress-equivalent-strain input data from subroutine 
ELAS until the occurrence of initial yielding. The stress- 
strain solution then proceeds along the yield surface as 
determined from the stress-strain properties in subroutine 
YIELD. At each increment during yielding the stress shift 
(difference between new yield stress and stress predicted 
from elastic analysis) from the original input data is 
calculated. Elastic load reversal is signaled when the 
input stress is less than the yield stress from the previous 
increment. During elastic unloading the stresses are 
translated from the original elastic-analysis solution by 
the amount of the calculated stress shift. Reverse yielding 
occurs when the stress reaches the reverse yield surface as 
determined from the hardening model incorporated in 
subroutine YIELD. Again the solution follows the yield 
surface until another load reversal is indicated when the 
stress, based on the shifted elastic solution, is less than 
the yield stress. The elastic response during load reversal 
is obtained by translating the original elastic solution 
according to the new stress shift calculated during 
reversed yielding. The stress-strain response for 
subsequent cycles is computed by repeating this 
procedure of identifying load reversals, tracking reverse 
yield surfaces, and translating the original elastic solution 
during elastic loading and unloading. Creep 
computations are performed for increments involving 
dweii times by using the creep equation and strain 
hardening rule incorporated in subroutine CREEP. 
Depending on the nature of the problem, the creep effects 
are determined on the basis of one of the three options 
provided in the subroutine. 

The code automatically avoids the Neuber-type strain 
redistribution corrections for thermal loading problems 
where there are no applied mechanical loads. Without 
applied loads, the Neuber method would be inapplicable 
since equation (4) w r ould have zero net stresses and strains 
in the denominator. Provision is also made for the user to 
circumvent the corrections for other situations where 
they should not be taken into account. These situations 
include locally strain controlled problems and, as shown 
in the next section, problems where the total strain input 
is based on strain measurements rather than elastic finite- 
element analyses. 

Since the stable cyclic stress-strain curve is a function 
of the plastic-strain range, it is necessary to iterate 
between the maximum plastic strain used to specify the 
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stress-strain curve and the analytical solution. In both the 
MARC code and the original ANSYMP code the 
iteration required rerunning the problem until the 
assumed and calculated plastic strains were in reasonable 
agreement. In the improved version of the ANSYMP 
program, the iteration was automated so that the analysis 
was repeated with the previously calculated maximum 
plastic strain used to define the stress-strain curve. This 
iterative process is continued until the specified and 
calculated maximum plastic strains agree within 10 
percent, usually within three iterations. 

The computer program was verified by conducting 
simplified analyses for the benchmark notch problem and 
comparing the results to those from MARC nonlinear 
analyses. The geometry of the benchmark notch 
specimen is illustrated in figure 3. This specimen was 
tested under isothermal conditions as part of a program 
to provide controlled strain data for constitutive model 
verification (ref. 5). A MARC analysis of this problem 
using kinematic hardening demonstrated excellent 
agreement with experimental data in reference 8. Two 
variations of this problem were analyzed, one with the 
initial loading carried out to a plastic strain of 0.4 percent 
and the second with the specimen loaded to a maximum 
plastic strain of 0.6 percent with a small amount of strain 
ratcheting. A kinematic hardening model was used with 
cyclic stress-strain data for Inconel 718 alloy obtained 
from reference 5. Nonlinear and elastic MARC analyses 
of this problem were performed by using approximately 
600 triangular-plane strain elements to model a quarter 
segment bounded by planes of symmetry (fig. 4). The 
MARC solutions shown for the benchmark notch 
specimen were computed at the closest Gaussian 
integration point to the root of the notch. Excellent 
agreement between the MARC solution and the measured 
strains at the notch root was demonstrated in the work of 
reference 8. 


71. 12 (radius) -v 



Figure 3, — Benchmark notch specimen [K, = 1.9). (Dimensions in mm.) 



Figure 4.— Benchmark specimen finite-element model. (Dimensions in 
mm.) 


Discussion of Analytical Results 

The results of the improved simplified inelastic 
analyses of the benchmark notch problems are discussed 
in this section. Comparisons are made with MARC 
inelastic solutions. Stress-strain cycles used for 
comparison purposes are discussed in terms of effective 
stresses and equivalent total strains by using the von 
Mises yield criterion with signs computed from the 
assigned signs in the original elastic solution. The entire 
discussion is based on the critical location at the notch 
root. 

The benchmark notch test was conducted by 
mechanical load cycling at a constant temperature of 649 
°C. A mechanically loaded structure, especially where the 
peak strain occurs at a stress raiser, is most likely to 
violate the basic assumption of the simplified approach 
that strain redistribution is prevented by containment of 
the local plastic region by the surrounding elastic 
material. Reference 1 reports that the total strain range 
from the MARC elastic-plastic analysis was 20 percent 
greater than that obtained from the MARC elastic 
analysis. This foreshortening of the elastic-strain range 
caused the simplified procedure to truncate the stress- 
strain hysteresis loop (fig. 5(a)). When the input total 
strain history was based on optical strain measurements 
at the notch root, the agreement was good between the 
simplified and MARC elastic-plastic stress-strain 
hysteresis loops, as demonstrated in figure 5(b). Both the 
ANSYMP and MARC elastic-plastic analyses gave stable 
stress-strain hysteresis loops for the second cycle. Note 
that when the simplified analysis uses strain 
measurements rather than elastic finite-element 
solutions, strain redistribution corrections are 
unnecessary and should not be applied. 

In figure 5(c) similar comparisons are shown by using 
the simplified procedure corrected tor strain 
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(a) Elastic solution as input, uncorrected for strain redistribution. 

(b) Strain measurements as input, uncorrected for strain redistribution. 

(c) Elastic solution as input, corrected for strain redistribution on loading part of cycle. 

(d) Elastic solution as input, corrected for strain redistribution on loading and unloading parts of cycle. 

Figure 5. — Comparison of ANSYMP and MARC results for benchmark problem (0.4 percent maximum plastic strain). 


redistribution on the loading part of the cycle and the 
MARC elastic solution as input. Neuber corrections for 
residual stresses and strain redistribution on the 
unloading part of the cycle were no! implemented for this 
case. The truncation of the stress-strain hysteresis loop 
(fig. 5(a)) was eliminated, and the ANSYMP solution 
agreed well with the MARC solution. This demonstrates 
the significant improvement in accuracy that can be 
attained by applying the Neuber correction in the plastic 
region, even without the residual stress correction during 
unloading. 


When both the loading and unloading strain 
redistribution corrections were applied, the agreement 
between the predicted and MARC elastic-plastic results 
was even better (fig. 5(d)). 

The original benchmark problem analyzed in figure 5 
was a completely closed cycle and the maximum plastic 
strain was about 0.4 percent. To exercise the improved 
simplified procedure on an even more severe case, the 
mechanical loading was increased so that the plastic 
strain reached approximately 0.6 percent, and strain 
ratcheting was induced. Again the agreement between the 
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liiune 6. — Comparison of ANSYMP and MARC results for benchmark problem (0.6 percent maximum plastic strain). 


ANSYMP and MARC inelastic solutions was excellent 
(fig. 6). Note that the simplified procedure was able to 
capture the strain ratcheting on the second cycle. 

The slight discrepancies between the the ANSYMP and 
MARC hysteresis loops in figures 5(d) and 6 were due 
primarily to the iteration process built into ANSYMP, 
which resulted in the use of a cyclic stress-strain curve 
slightly different from that used in the MARC elastic- 
plastic analysis. The ANSYMP analyses of the 
benchmark notch problem used 0.3 percent of the central 
processor unit time required by the MARC nonlinear 
analyses; this is without considering that the latter had to 
be run several times to obtain a cyclic stress-strain curve 
compatible with the calculated maximum plastic strain. 

Summary of Results 

An improved simplified analysis procedure was 
developed for calculating the stress-strain history at the 
critical location of a thermomechanically cycled 
structure. This improved procedure incorporated 
Neuber-type corrections to account for strain 


redistribution and residual stresses due to plastic-strain 
reversals. The following results were obtained from the 
evaluation of the method: 

1. The predicted stress-strain response based on elastic 
finite-element solutions with the Neuber-type corrections 
showed excellent agreement with elastic-plastic finite- 
element solutions using the MARC program. 

2. The predicted stress-strain response using only the 
corrections for plastic yielding on the loading part of the 
cycle showed good agreement with the MARC elastic- 
plastic solution. The corrections for the unloading part of 
the cycle were of secondary importance. 

3. The strain redistribution corrections should not be 
applied when the input total strain history is based on 
local strain measurements rather than elastic finite- 
element analyses. It is also known from previous 
evaluations of the simplified procedure that these 
corrections should not be applied in thermal loading 
problems. For cases where a structure is subjected to a 
combination of thermal and mechanical loads, a method 
will have to be developed to partition the thermal and 
mechanical stresses so that the strain redistribution 
corrections are applied only to the mechanical stresses. 
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4. Stress-strain hysteresis loops were computed at the 
critical location of the structure using 0.3 percent of the 
central processor unit time required for elastic-plastic 
finite-element analyses. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, January 8, 1985 
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